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In this paper we study the accretion of dark energy onto a black hole in the cases that dark 
energy is equipped with a positive cosmological constant and when the space-time is described by a 
Schwarzschild-de Sitter metric. It is shown that, if confronted with current observational data, the 
results derived when no cosmological constant is present are once again obtained in both cases. 

PACS numbers: 98.80.Qc, 03.65.Fd. 



I. INTRODUCTION 

In the last few years a large number of models have 
been proposed to describe the current accelerated expan- 
sion of our Universe [l| . The consideration of these new 
models entails the appearance of new and surprising phe- 
nomena, such as those named big rip Q , big trip [f| or big 
hole [g, ID; la] ■ One of the most popular models describing 
the accelerated expansion state is the quintessence model 
. In this scenario one considers a time dependent scalar 
field which can be interpreted as an homogeneous and 
isotropic fluid with an equation of state, p = wp (where 
units, c = G = 1, are used), in which p and p are the 
pressure and energy density, respectively, and w is a con- 
stant which current observations make to look close to 
— 1, but with a given bias toward the phantom regime 
H, w < —1. The latter case corresponds to a "super- 
accelerated" expansion of the Universe, until it finally 
reaches a big rip singularity {ij. If present wormholcs 
accrete that phantom energy, then their sizes increase so 
rapidly that can eventually engulf the universe itself, so 
allowing for cosmic time traveling [f| . On the other hand, 
present black holes accreting dark energy with w > — 1 
could in principle eventually engulf the universe in a pro- 
cess that has been dubbed big hole [6]. In this case, 
nevertheless, current observational data prevent such a 
process from occurring |6j. 

However, a question which still remains unanswered 
is what happens in the case in which the cosmic expan- 
sion is caused by a certain dark stuff when a positive 
cosmological constant is involved at. In that case, could 
the accretion of dark energy onto black holes be so fast 
that the Universe would eventually undergo a big hole 
phenomenon?. In order to answer that question, two dif- 
ferent ways can be followed. The first one is developed 
in sec. II, where a cosmological constant is introduced 
in the Friedmann equations, taking for the mass rate of 



the black hole the expression obtained from a non-static 
Schwarzschild metric [f|. In that case, if we consider 
an homogeneous and isotropic model with a black hole 
described by a Schwarzschild metric, no big hole phe- 
nomenon would happen Q . In the second way Friedmann 
equations with no cosmological constant are considered, 
using the Schwarzschild-de Sitter metric to describe the 
space-time. This allows us to study in sec. Ill the accre- 
tion procedure in a cosmic scenario, where the universe is 
no longer asymptotically flat but asymptotically de Sit- 
ter. Finally, in sec. IV, some conclusions are drawn and 
further comments added. 



II. ACCRETION OF DARK ENERGY WITH A. 

Let us consider first a dark energy model described 
by a perfect fluid with an equation of state given by 
p = wp, where w is a constant parameter. Integrat- 
ing the expression for the cosmic energy conservation, 
dp = — 3(p + p)da/a, where a(t) is the scale factor, we 
obtain 



p = p (a(t)/a )- 3 ( 1+w \ 



(1) 



In that case, the Friedmann equation for a flat universe 
with a cosmological constant reads, 



\ + Ca- 3< - 1+w \ 



(2) 



in which the overhead dot denotes the derivative with 
respect to time, A = A/3 > and C = 8np /(3aQ 3< " 1+w ^), 
with clq = a(to), being to the time at the current epoch. 
Eq. ([2]) can be integrated to give the time evolution of 
the scale factor @, Q • We get 
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where 



D = 



(4) 



from which it can be easily checked that, < D < 1. 

As is well known, the general expression for the black 
hole mass rate in the case where the black hole accretes 
dark energy in a non-static Schwarzschild space, is given 
for an asymptotic observer by [||, 



M = 4irAM 2 (p + p), 



(5) 



where A is a constant. This rate equation can be inte- 
grated in any homogeneous and isotropic model, giving 



M(t) 



Mn 



(6) 



1 + AM [H(t) - H ] 
which, taking into account Eq. Q, can be re-expressed 



M{t) = M 



2A 1 / 2 DA/nA e 3 ( 1 +«') Al/2 ( t - t ») 
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1 - D 



e 3(l+ lU )Ai/2( t _ to ) _ jj 



If the Universe is filled with a phantom fluid, for which 
w < —1, the scale factor given by Eq. (|3|) diverges at 
the so-called big rip singularity [J], which takes place at 
a finite time tbr in the future given by 



tbr — U) + 



|Ln.D| 



3(M - ljA 1 ^- 



(8) 



From Eq. ([7]), it can be seen that the black hole mass 
can be written as, M = M [l — F(£)] -1 , where F(t) is 
a decreasing function of time and F(t ) — 0, i. e., F(t) 
is a negative function. Then, the black hole mass will 
decrease progressively until it vanishes at the Big Rip, 
for any Mo. Therefore, any black hole in a phantom 
dominated universe tends to disappear at the big rip sin- 
gularity, a result that was already derived by Babichcv 
et al. [10] for the static case with A = 0. 

On the other hand, if a fluid with w > — 1 is consid- 
ered, then the scale factor of the Universe would steadily 
increase forever (see Eqs. Q and ([?])). In that case, 
the function F(t) becomes an increasing function of time 
with a finite limit given by, 7 = 2A 1 / 2 M IL4/(1 - D), 
when t — ► oo. Then, the Universe will suffer a big hole 
phenomenon if 7 > 1. Expliciting the expression for D 
leads to the condition 



M A [(A + 87rp /3) 1/2 -A 1/2 



> 1. 



(9) 



However, from the WMAP data, it follows that Ho ~ 
10 _26 (meters) _1 , and making use of Eq. @, it can be 
estimated that (A + Sirpo/3) 1 / 2 ~ ICU 26 , and generally, 
7 ~ Mp (lfT 26 - A 1 / 2 ) < M 1CT 26 < 1, as it is shown 
in [1, Therefore, there could exist no black holes in 
the Universe with a initial mass so big as to produce 
a big hole phenomenon. The conclusion can be more- 
over drawn that for a black hole to produce a big hole 
phenomenon in a universe with a positive cosmological 
constant it is necessary that originally it be even bigger 
than those that could produce a big hole in a universe 
without cosmological constant. 



III. ACCRETION ONTO SCHWARZSCHILD-DE 
SITTER. 

In the Babichev et al. model |10], and its general- 
ization for the non-static metric [6[, the black hole was 
described by means of the Schwarzschild metric (its non- 
static generalization), immersed in dark energy depicted 
as a perfect fluid. In the previous section we have inserted 
a positive cosmological constant in the Friedmann equa- 
tions. We had therefore a universe filled with an "effec- 
tive fluid" containing both dark energy and a cosmolog- 
ical constant, where a static or non-static Schwarzschild 
black hole was immersed. 

In the present section it will be instead assumed that 
the stuff which produces the expansion of the Universe 
is a dark energy fluid alone, and the contribution of 
the cosmological constant is taken into account through 
the space-time metric itself; that is, we consider now 
a Schwarzschild-de Sitter metric to describe that space- 
time, 



ds 2 



where, 
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(10) 



(11) 



The possible horizons of this space would be the zeros of 
the A function. One can see that if 3MA 1 / 2 > 1 there 
is only a real solution which is negative; if 3A/A 1 / 2 = 1, 
there is one negative real solution and two equal positive 
real solutions; and if 3M A 1 / 2 < 1, then there is one nega- 
tive real solution and two different positive real solutions. 
In the latter case, one can provide these two positive real 
solutions with the physical meaning of two horizons [ll[ • 
So, in that case, there is a black hole horizon r^h and a 
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cosmological horizon r c , defined by 
2 



A 1 / 2 
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and 



A 1 / 2 



cos 



arccos(-3MA 1 / 2 ) 



(12) 



(13) 



If the black hole mass M increases, the size of the 
black hole horizon would also increase whereas that of 
the cosmological horizon would decrease, until a limiting 
case in which M takes such a value that 3MA 1 / 2 = 1, 
where the two horizons reduce to just one horizon, so 
disappearing the physical space between them. 

Now, let us assume that the space is filled by dark en- 
ergy which is described by a perfect fluid with an energy- 
momentum tensor given by, 



Tya, = (p + 9)11^ - pg^, 



(14) 



where p and p are the pressure and energy density, respec- 
tively, and is the four- velocity, such that u IL u ^ = — 1. 

We consider, such as it was made in Refs. [l(| and @, 
the time component of the conservation law of the energy- 
momentum tensor, Tq. v — 0, and derive the following 
equation 



In order to integrate Eqs. (jl5)) and (|T7|) over the ra- 
dial coordinate one must carefully choose what are the 
integration limits. In principle, we would like to consider 
an observer placed on the region rj,h < r < r c . Since 
there are two horizons, one can integrate Eqs. (|15[) and 
(fl~7|) between two different pairs of limits: i) from the 
cosmological horizon up to the observer, or ii) from the 
black hole horizon up to the observer. In the first case 
one obtains, 



ur 2 [p + p) [u 2 + ,, 



1/2 



exp 



drf(r,t) 



and 



Ci(t) 
(19) 



r it exp 



dp 

p + p 



exp 



drg(r,t) = B 1 (t), 



(20) 

with C\{t) — —r 2 (p + p) \ r =r a and Bi{t) = — r 2 , where 
we have considered u — ► — 1 when r — > r c . On the other 
hand, if the lower integration limit is placed at the black 
hole horizon, then we have 



ur (p + p) I u + 



1/2 



exp 



drf(r,t) ) =C 2 (t), 
(21) 



where 



d r T5 + -T£ + f(r,t) = Q, 



f(r,t)=d T° + — (T;-T°) 



(15) 



(16) 



On the other hand, the projection onto the four-velocity 
of the energy-momentum tensor conservation law gives, 



d r p d r u 2 , . 
-j L + — + ~+9(r i t) = 0, 
p + p u r 



with u = u r and 



9{r,i) 



u°dop dou 



u (p + P) u 



(17) 



(18) 



r u exp 



where (^(i) 




' bh 



exp 



(P + P) 



drg(r,t) = B 2 (t), 



and B 2 (t) 



(22) 

-r 2 
'bh' 



with a similar consideration on u and r. Obviously, Eqs. 
(fT9|) and ([2T|) must be equivalent, as it should also hap- 
pen with Eqs. (|2"01) and (f2"2"| . In what follows it will be 
seen that this is actually the case, because the same re- 
sult is obtained by evaluating Eq. (fT9)) (or (20)) at the 
black hole horizon as Eq. (|2"Tj) (or (22)) at the cosmologi- 
cal horizon. It will be also shown that the integral of the 
functions /(r, t) and g(r, t) must be finite as it should be 
expected. 

Taking into account Eqs. ([T9| and (|20| it can be ob- 
tained, 



(p + p)[vr + -*) exp - 



1/2 



dp 

Pr c P + P, 



exp / dr [f(r, t) - g(r, t)} = (p + p) \ r=rc , 



(23) 



and from Eqs. (J3TJ) and (|2^|) . 
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{p + p)\u z + 



1/2 



exp 



dp 



exp I / dr[/(r, t) - g(r, t)] = (p + p) \ r=rbh ■ 



(24) 



r 



Due to the spherical symmetry, the rate of black hole 
mass can be expressed as M = Attv^Tq, or using Eqs. 
PPf and ((231) or Eqs. (HH and QZQ, as 



for an observer at the black hole horizon, which has the 
same right-hand-side as in Eq. (j29|) . This can be easily 
integrated in a quintessence model, where 



M = Airr 2 (p + p) | r=rc exp 



M = Anr 2 h (p + p) \ r=rbh exp 



drf(r,t) 



(25) 



drf(r,t) . (26) 



Thus, evaluating Eq. (|2"5|) at r c and Eq. (f2"o) at r&h, one 
easily obtains the black hole mass rate for observers at 
the cosmological horizon and at the black hole horizon, 
respectively, which turn out to be 



M = Airr 2 c (p + p) | r=rc 



(27) 



and, 



p = p [a(t)/a r 3{1+w) 



Po 



l + -C{l+w)(t-t ) 



With C now defined by C = (87rp /3) 



1/2 



(31) 



Using Eqs. (12) and (13), we can integrate the left- 
hand-side of Eqs. (29) and (30) to give 



1 [A ( 3 + 2z 3 + 2z 



2 V 2 V\/TTI y/T+zZ 
where z is given by, 



4irp (l + w)(t-t ) 

l + %(l + w)C(t-to)' 
(32) 



M = 4irr 2 bh (p + p)\ r 



(28) 



Therefore, in order to obtain the black hole mass func- 
tion in an isotropic and homogeneous model one must 
integrate the following equation 



dM , , 

— = 4tt (p + p) dt, 



(29) 



z = cos {2/3 arccos(~3MA 1/2 ) + 4tt | , (33) 
in the first case, or by 



cos 2/3arccos(— 3MA 1 / 2 ) 



(34) 



for an observer at the cosmological horizon, and 

dM A . 

-2- = 4tt (p + p) dt 

r bh 



(30) 



in the second one. Expressing the mass of the black hole 
in terms of z by the inverse of the previous expressions, 
i.e. using 



x(z) = 3VAM(z) 



- arccos z — 4ir } = 
2 



cos - arccos z 
2 



(1 - 2z) 



l + z 



(35) 



r 



one can note that the two different definitions of z, given 
by Eqs. (|33[) and (|3~4"|) , which imply two different equa- 
tions expressed through Eq. ([32|) . actually produce the 
same mass function, i. e., the black hole mass function 
for an observer at the cosmological horizon is the same 
that for an observer at the black hole horizon. There- 
fore, one could assume that this is a well-defined problem 
and that there is a unique black hole mass function inde- 
pendent of the placement of the observer. On the other 
hand, because z arises from an arccosine function in Eq. 



(|35| . it might at first sight range from z = — 1 to z = 1, 
but actually it has to take only on values in the interval 
— 1 < z < i to keep on non- negative values for the mass 



function (|35[) . In fact, x(z) is an increasing function from 
z = — ltoz = — h, where it reaches the value x = 1. It is 
worth noticing that it is precisely the latter value the one 
at which both horizons coincide. Thereafter, x keeps be- 
ing a decreasing function all the way until z = 5, where 
x vanishes, i.e, when the black hole has evaporated away. 

As we have already mentioned above, the condition 



5 



that allows us to have a physical space between the two 
horizons is 3MA 1 / 2 < 1. Therefore, this condition may 
be considered to be also fulfilled at the current time where 
the black hole has a mass Mo, i.e., when xq < 1 , where 
x = 3A/A 1 / 2 . In this way,at least in principle, a Taylor's 
expansion could be a good treatment for Eq. (j32|) . In 
that case, one has 



where the function A(t) is given by, 



x 



4 

3^0 



with 



F(t) 



3/2(1 + w)C(t- t ) 
1 + 3/2(1 + w)C{t - t ) ' 



(36) 



(37) 



Note that F(t) would be an ever increasing function with 
time if 1 + w > 0, and a decreasing function with time 
if 1 + w < 0. It follows then from Eqs. (J^ZJ) and f2"g]) 
that black holes would increase their size if they accrete 
dark energy with w > — 1, and would shrink if they ac- 
crete phantom energy with w < —1. Therefore, in the 
phantom case x will take always on smaller and smaller 
values as the time grows up, making Eq. (|36p more and 
more accurate as the the black hole is evaporating away 
to disappear at the big rip time 



tbr — to 



1 



(\w\ - 1)^/671-/90 



(38) 




h, t 



FIG. 1: Time evolution for the mass of a black hole within 
a Schwarzschild-de Sitter universe filled with a phantom fluid 
(w = -1.1), as evaluated by x(t) = 3A 1/2 M(t), for different 
initial masses. From top to bottom: (a) xq — 2.1 10~ 3 , (b) 
xo = 1.4 10~ 3 and (c) xo = 710" 4 , in units Hq 1 , i.e., Ho = 
= 1 and H < A = 1.1. 



The previous line of reasoning is based on an approxi- 
mate method. One could also try to perform an analyti- 
cal study in closed form by re-expressing Eq. ([32]) as, 



A(t) = A a - 



cg(t - t ) 
l + Cl (t-t ) : 



with 



3 + 2z 



^0 = 

c = 87rp o y ^ (1 + w), 
ci = ^6%po (1+ w), 



(40) 



(41) 

(42) 
(43) 



and z = cos (| (arccos(— 3AfoA) + 4tt)), M being the 
initial mass of the black hole. If we take io « 1, which is 
a sufficiently good approximation for any practical initial 
mass of the black hole, then zo ~ \ (the value zo = — 1, 
which corresponds to an initial black hole mass very close 
to cero, is neglected in order to preserve a well-defined 
meaning for Ao). 




FIG. 2: Time evolution for the mass of the black hole within a 
Schwarzschild-de Sitter universe filled with a dark energy fluid 
(w = -0.9) as evaluated by x(t) = 3A 1/2 M(i), for different 
values of the initial mass: xo = 0.88, xo = 0.70 and xo = 0.58. 

Squaring Eq. ((39|) we can attain the following two solu- 
tions, 



z± 



A 2 {t) - 12 ± \A(t)\y/A 2 (t)-8 



8 



(44) 



3 + 2z = A(t)VT+7„ 



(39) 



These two values, z±, inserted in Eq. ([55]) . imply two 
solutions for the time evolution of the black hole mass. 
The solution z + yields a negative value for the initial 
mass of the black hole, both in the phantom and in the 
dark energy regime, so it ought to be disregarded. The 
mass evolution for the value z_ is depicted in Figs. 1-3. 
In the phantom case, the time left until the big rip singu- 
larity turns out to be, th r — tn — -n- — , 1 .„ Then, from 

Eqs. (|4"0"]) and (|4"T]) . it can be checked that as the uni- 
verse approaches the big rip, the function A(t) tends to 
infinite, i.e., A(t — > tbr) ~ > irrespective of the initial 
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mass of the black hole. Therefore, because in that limit, 
z_ — ► — 1 and x(z_) —> 0, any black hole eventually ends 
up completely evaporated at the time in which the big 
rip singularity takes place. That is the same conclusion 
as in Refs. [UJ and [8[- It follows that in the phantom 
case it does not make any difference between consider- 
ing the black hole in an asymptotically flat scenario or 
in a cosmological scenario (see Fig. 1 which corresponds 
to different values of the initial mass of the black hole). 
On the other hand, Fig. 2, gives the case for dark en- 
ergy, w > — 1, where the mass of the black hole grows 
up, and therefore the size of its event horizon also in- 
creases whereas the cosmological horizon decreases until 
it eventually reaches the cosmological one, at x = 1. 




FIG. 3: Time evolution for both the black hole horizon, rbh, 
and the cosmological horizon, r c , Eqs. (|12|l and (113[) . for 
different values of the parameter w. From left to right: wi — 
-0.88, W2 = -0.92 and w 3 = -0.95, with A = 3.5 and the 
same units as those used in Fig. 1. 



The time at which both the cosmological and the black 
hole horizons coincide depends on the value of the param- 
eter w. 



t,n — tr 



where, C = Ag 2V2 • constant. This 

case is depicted in Fig. 3. As the value of w approaches 
— I, the time at which the universe is engulfed by the 
black hole grows up. Therefore, if we consider a cosmo- 
logical scenario, where the black hole space is asymp- 
totically de Sitter a big hole phenomenon would take 
place, which is so derived not only for the black hole 
horizon growth but also because the cosmological hori- 
zon decreases. As it should be expected in the limit of 
w = —I, the size of the black hole remains constant in 
time, and the big hole phenomenon does not take place. 



C 



(45) 



IV. CONCLUSIONS. 

In this work we have analyzed the accretion phe- 
nomenon of dark energy onto non-static black holes in 
a consistent cosmic scenario. As a first approach, we 
have considered a non-static Schwarzschild black hole 
onto which dark energy is accreted in the presence of a 
positive cosmological constant. As it was already shown 
in Ref. Q, upon accreating phantom energy, any black 
hole decreases its mass and tends to disappear at the 
big rip singularity. Coversely, if the parameter w in the 
equation of state is greater than —I, the black holes in- 
crease their sizes in a moderate, safe way since they can 
never reach a size large enough to be able to engulf the 
universe. 

Then, in order to consider a more realistic cosmic sce- 
nario for dark energy accretion, we have taken a non- 
static Schwarzschild-de Sitter universe, in which the mass 
is an arbitrary function of time. In this cosmological con- 
sistent framework, the accretion of phantom energy onto 
black holes makes the size of the hole to steadily decrease 
to finally vanish at the big rip time, such as it is shown in 
Fig. 1. We would like to remark that the same result was 
already obtained in Ref. [Io| , where this phenomenon 
was studied using a static Schwarzschild metric. The rea- 
son for this agreement is that the accretion rate obtained 
in the mentioned work may only be valid at small accre- 
tion rates, a regime that is always satisfied in case of a 
phantom fluid. 

For w > —I in the mentioned cosmologically consis- 
tent scenario, accreating black holes would increase their 
masses. The growth in the black hole mass would pro- 
duce an increase in the black hole horizon size, and a 
decrease in the cosmological horizon size until the two 
horizon join up into just one, when M(t) = 1/(3A 1 / 2 ). 
Therefore, a big hole phenomenon would take place in 
this case because the whole universe, delimited by the 
cosmological horizon, is within the black hole horizon. 

It is worth noticing that in this paper only a classical 
treatment has been performed. It would be expected that 
quantum effects would probably smooth these phenom- 
ena, leading to a non- vanishing black hole mass at the big 
rip in the phantom case, |12J, and, in the case w > —1, to 
a deceleration in the black hole horizon growth and in the 
cosmological horizon decrease, due to particle creation at 
the two horizons induced by the Unruh effect. 

Finally, we would like to point out that in a recent 
work Gao et al. [IH declared that the mass of a black 
hole cannot decrease due to phantom energy accretion, 
in contradiction with the results achieved in the present 
work. We consider that their result is wrong because 
they introduced a premise in which their previously de- 
sired result is actually contained, making circular their 
whole argument. In particular, these authors choose a 
particular function of the black hole mass proportional 
to the scale factor at the beginning of section III and at 
the end of section IV. Obviously, if one assumes a black 
hole mass function proportional to the scale factor, the 
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resulting black hole mass would only depend on the expo- 
nent of the scale factor entering the mentioned function 
and on the time derivative of the scale factor (which is 
positive for all expanding universes). On the other hand, 
their initial assumption is totally contradictory with all 
studies on dark energy accretion onto black holes in the 
literature for which the time rate of the black hole func- 
tion is recovered at least qualitatively at some limit, such 
as it occurs in the present work. 
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